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Abstract 

We numerically study the distribution function of the conductance (transmission) in the one- 
dimensional tight-binding Anderson and periodic-on-average superlattice models in the region of 
fluctuation states where single parameter scaling is not valid. We show that the scaling properties 
of the distribution function depend upon the relation between the system's length L and the length 
l s determined by the integral density of states. For long enough systems, L>/ s , the distribution 
can still be described within a new scaling approach based upon the ratio of the localization length 
l\ oc and l s . In an intermediate interval of the system's length L, l[ oc <C L <C l s , the variance of the 
Lyapunov exponent does not follow the predictions of the central limit theorem and this scaling 
becomes invalid. 

PACS numbers: 72.15.Rn,42.25.Bs,41.20.Jb 
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I. INTRODUCTION 



Coherent transport properties of disordered systems have been a subject of active research 
for the last thirty years, but complete understanding of this phenomenon even for one- 
dimensional models is still absent. Even though the scaling theory, put forward in the 
pioneering work of Ref. U, created a successful conceptual framework for discussing the 
phenomenon of localization, the theoretical foundation of the scaling hypothesis itself has not 
yet been completely understood. One of the principal difficulties that the scaling theory of 
localization had to deal with from the very beginning was an absence of self-averaging of the 
main transport coefficients: conductance, g, or transmission, T. Therefore, even the nature 
of the scaling parameter remained unclear until it was realized that the scaling hypothesis has 
to be applied to the entire distribution function of the conductance or transmittance.^ 4 ^ 

For one- dimensional systems Anderson, et al.— suggested that the most suitable quantity 
for dealing with the statistical description of conductance is the Lyapunov exponent (LE), 
which can be defined for systems with finite length L as 



The name 'Lyapunov exponent' alludes to the fact that the quantity defined by Eq. (0) have 
the same statistical properties as the "real" Lyapunov exponent, i.e. the exponential growth 
rate, {1/L) log of the norm of the wave function, ip. An important property of LE is 
that it satisfies a multiplicative central limit theorem^ and approaches a non-random limit, 
7, when the size of the system, L, tends to infinity. The localization length, Z/ oc , of a state 
with energy E in the infinite system is related to 7 as l[ oc = 7 -1 . At finite L, 7 is a random 
quantity with mean value equal to 7, (7) = 7. The distribution of LE is the main object 
of research in the field of one-dimensional localization. The hypothesis of single parameter 
scaling (SPS) in this context means that the distribution function can be parameterized by a 
single parameter, 7 itself. As a result, it is expected that all moments of the distribution can 
be expressed in terms of the first moment, (7), in a universal way. For the second moment 
(variance), a 2 , such a relationship, as it was first conjectured by Anderson et al.,— can be 
presented in the form 

° 2 - 1 < 2 > 

The entire distribution function of LE for systems with finite lengths was also derived 




(1) 
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by several authors in the limit of infinitesimally weak local scattering for several models.— 
For finite L, this function was found to be non-Gaussian, but nevertheless, it depended upon 
a single parameter - the localization length. 

Thus, in the situations when SPS holds the problem of the conduct ance/transmision dis- 
tribution function can be considered as settled. There are spectral regions, however, where 
SPS fails even for locally weak disorder. These are, first of all, the regions of fluctuation 
states, which arise outside of the initial spectrum because of disorder. This result was first 
obtained numerically in Ref. Q for a periodic-on-average system and was confirmed by an 
exact analytical solution of the Lloyd model (the Anderson model with the Cauchy distri- 
bution of the site energies).—' 12 Similar results were obtained numerically for the Anderson 
model with the bo x 11 ! 12 and dichotomic 13 distributions of the site energies, and analytically 
for a continuous model with white noise Gaussian potential.— The analytical calculations of 
Refs. Illll2 revealed that the criterion for the validity of SPS can be presented in the form 
hoc > h, where l s is a new scale introduced in Refs. 1 1 ill 1 2t For the Lloyd model this scale 
is defined in terms of the imaginary part of the Lyapunov exponent, which, according to 
Thouless^ is proportional to the integral density of states. Therefore, l s can be presented 
in the form 

h = sin (irN(E)Y (3) 
where N(E) is the integral density of states between the genuine boundary of the spectrum 

and the energy E normalized by the total number of states in the system. The definition of 
l s in this form can be easily generalized to other models as well, and it was shown numer- 
ically that the SPS criterion based upon l s works for such models as the Anderson model 
with bo x 11 ' 12 and dichotomioi 3 - distributions of site energies, a Kronig-Penney-like model 
with a periodic-on-average distribution of barrier widths , 11112 and a model of a scalar wave 
propagating in an one-dimensional absorbing disordered medium.— The case of periodic- 
on-average models involves a system with multiple bands, and in this case N(E) must be 
understood as the integral density of states between a genuine boundary of the band (if the 
latter exist) and the energy E normalized by the total number of states in the band. More 
detailed discussion of this case can be found in Ref. 12 . In a recent paper Ref. Q, it was 
shown how this criterion can be applied to the zero energy states of the Anderson model 
with a diagonal disorder, where the violation of SPS was observed in Ref. Q. 

The criterion based on l s replaces an original criterion put forward by Anderson, et al.— 



that suggested that SPS exists if the stationary distribution of the phases of the reflection 
and transmission coefficients is uniform, and the phase relaxes to this distribution over a 
length, which is much smaller than the localization length. By using the hypothesis of the 
phase randomization Eq.|2]was re-derived by many authors for a variety of different models. 7 
The phase randomization was proven rigorously in some one-dimensional ^ 19 ' 20 and quasi-one- 
dimensiona l 21 ' 22 models, but only for certain parts of the spectrum of the respective systems. 
At the same time, it was found that, for instance, in the Anderson model with a diagonal 
disorder the stationary distribution of the phase is not uniform for all values of energy E, for 
which cos -1 (E/2) is a rational fraction of ir (it is assumed that in the non-random case all 
site energies in the Anderson model are set to zero, and the interaction parameter is chosen 
to be equal to unity). The strongest deviation of the phase distribution from the uniform 
one takes place in the vicinities of E = and the initial band boundaries E = ±2. While it 
was found that an absence of the phase randomization in both of these cases is accompanied 
by the violation of SPS , 10 ' 11 ' 12 ' 18 the reference to the phase randomization as a criterion for 
SPS does not seem to be satisfactory. Indeed, the initial idea of the phase randomization 
length,— used to introduce the criterion for SPS, does not actually describe the way the 
distribution of phase becomes non-uniform. The absence of the phase randomization does 
not mean that the relaxation length of the distribution of phase becomes too large and 
exceeds the localization length. What it means is that the stationary distribution of phase, 
which can be reached over relatively short distances, is merely not uniform. Thus, the 
problem of a criterion for SPS is simply replaced by the problem of finding a criterion 
describing the transition between uniform and non-uniform stationary distributions of the 
phase. A solution for the latter problem suggested, for instance in Ref. [is, applies to only 



one particular model, and, actually involves different criteria for different spectral regions 
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121 was proven to work for the 



In contrast, the criterion based on l s introduced in Refs. 
entire spectrum of the variety of different models, and offers, therefore, a universal approach 
to the verification of SPS. 

The violation of SPS in the spectral region of fluctuation states rises a question about 
the properties of the probability distribution of LE in these regions. Recently, a significant 
progress in this direction was achieved in Refs Il4ll23l In the former paper, the first four 
moments of this distribution were found analytically for the Anderson model with a Gaussian 
white-noise potential. The authors of the latter paper used numerical simulations to develop 
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a macroscopic scaling approach to this problem, and one w 



rich could be readily applied to 



a wide variety of different systems. It was shown in Ref. 23, that not only second, but also 
the third moment of the distribution function of LE for the Anderson tight-binding model 
with diagonal disorder can be fully characterized by a scaling parameter k = h oc /h- 



The objective of the present paper is to present more fully and expand the results of 



Ref. 



23l Considering two quite different models of one- dimensional localization such as the 



Anderson tight-binding model with a diagonal disorder, and a model of a scalar wave prop- 
agating in a one-dimensional random superlattice, we demonstrate that one-dimensional 
disordered systems allows for a universal scaling description of the conductance (or trans- 
mission) distribution in the spectral regions of fluctuation states, where standard SPS does 
not work. In particular, we show that the scaling approach suggested in Ref. |23j describes 
not only the Gaussian bulk of the distribution function, but is also capable of describing the 
statistics of large deviations characterized by the third moment of the distribution. 

The results presented in this paper are also relevant to the problem of resonant tun- 
nelling through disordered potential barriers. For the first time, this problem was considered 
in the pioneering work by Lifshits and Kirpichenkov 2 ^ for quantum particles incident on a 
three-dimensional barrier, and was later studied in many subsequent papers (see reviews in 



Ref. 
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26T ). Mostly, these works were concerned with tunnelling through 3-dimensional bar- 



riers with the dimension in the propagating direction much smaller than in the perpendicular 
directions. Even though the resonant tunnelling is in many aspects a quasi-one-dimensional 
process,— the transport in the pure one-dimensional models significantly differs from the 
situation described above. First of all, in one- dimensional case all states are localized and 
transmission at any energy can be described as a resonant under-barrier tunnelling. There- 
fore, the difference between transport in the region of states from the initial allowed bands 
and the fluctuation states, is not as clear as in three-dimensional situations. Therefore, the 
problem of transport via fluctuation states was not considered as a separate problem in the 
area of one-dimensional localization until very recentl y 11112 Second, the main quantity of 
interest in the case of three dimensional barriers is the total transmittance across the entire 
area of the barrier, which is determined by the sum of individual transmissions through 
independent quasi-one-dimensional channels or filaments.— ^ This quantity approaches a 
non-random limit when the area of the barrier tends to infinity. In a pure one-dimensional 
case the self- averaging quantity is the Lyapunov exponent, which becomes non-random when 
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the length of the system becomes infinite. In a sense, the pure one-dimensional case is an 
opposite limit to the one considered for three-dimensional barriers. At the same time, so- 
lutions of the one-dimensional problem can be used to describe barriers whose lengths are 
larger than the typical localization length of individual channels. 

Another important application of the problem studied in this paper lies in the field 
of random lasing, which has become an area of active research.— It is anticipated that 
using localized modes of a strongly scattering disordered medium, one can obtain very low- 
threshold lasing. Disordered photonic crystals, which support fluctuation photon states in 
the band-gaps of the underlying periodic structures, can play an important role in achieving 
this objective . 28 ^ 9 The results presented in this paper will help to understand the unusual 
statistical properties of the lasing threshold and the nature of lasing modes in such structures. 



II. MODELS AND TECHNICAL DETAILS 



In this paper we study two models of one-dimensional Anderson localization: a classical 
Anderson tight-binding model with a diagonal disorder, and a scalar wave propagating in a 
one-dimensional random superlattice. The Anderson model is described by the equation of 
motion 

VWi + V'm-i + (U m - E)if) m = 0, (4) 
where random on-site energies U m are described by a uniform probability distribution: 



P(U„ 



The propagation of a scalar wave is described by a regular wave equation, 

d 2 ^ 




dx 2 



+ k 2 e(x)i) = 0, (5) 



with a piece-wise dielectric function, corresponding to a superlattice consisting of two types 
of layers with dielectric constants E\ and €2, respectively. The width of the layers of the first 
kind is kept constant and is equal to d\, while the width of the layers of the second type was 
chosen from a random distribution. In this paper, we report the results for (i) di uniformly 
distributed in the interval [d?) — 5, (d?) + 5 (uniform distribution), (ii) d<i taking one of two 
equally probable values — 5/^/3 and +5/^/3 (dichotomic distribution). 
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Both these models can be studied using the transfer matrix approach, in which the 
propagation of the excitation along the system is presented in the following form 

v m+1 = T m v m , (6) 

where two dimensional state vector, which presents the state of the system at the 

m-th site (or m-th interface between the layers) and T m is the transfer matrix describing 
the change of this state at one discreet step. For the Anderson model the state vector and 
the transfer matrix have the following forms respectively 

/ 



1pm 



(7) 



1 

For the second model the state vector can be defined as 

/ 



T m =\ E ~ Um ~ l \- (8) 



(9) 



A, 

ib' 

where tb m and tb' m are the values of the wave function, ip(x), and its derivative at the m-th 
interface between the layers. The transfer matrix in this case takes the form 

/ cos{k m d m ) (l/k m )sm(k m d m ) \ 
y-k m sm(k m d m ) cos(k m d m ) J 

where k m = ky/e^. The most important property of the transfer matrices is that the transfer 
matrix T M describing the evolution of the initial state vector across the M sites (slabs) is 
equal to the product of the one-step matrices 

M 

T M =n T «- ( n ) 

i 

Using the transfer matrices, we calculate the finite size LE, which for both models is defined 

as 

„ 1 ||T M tto|| , 10 v 
7 = yln^ — J]— , (12) 

where L characterizes the total length of the system. For the Anderson model, L — M if 
the distance between adjacent sites is chosen as a unit of length, and for the wave equation, 
L is a sum of the lengths of all slabs, and is a random quantity. 



We calculate LE iteratively using Eq. ()12j) starting with an arbitrary initial vector vq- 
The resultant vector is re-normalized after every ten iterations in order to avoid any loss of 
accuracy.— Since we are interested in statistics of finite size LE, we do not try to find its 
limiting value for L — > oo. Instead, we keep the size of the system fixed while calculating 
7 for different realizations of our systems. At the same time, since we are interested in 
asymptotic properties of the distribution, we consider only sufficiently long systems, for 
which L 3> lioc, where the localization length, li oc , is defined through the average value of 
LE as l ioc = (j)" 1 . 

Another quantity of interest in this work is the length l s , which is expressed in terms of 
the integral density of states N(E), Eq. (JHJ). For the Anderson model N(E) can be computed 
with the help of the node-counting theorem.— Starting with an arbitrary initial vector and 
the energy values, E < —2 — U, which are certainly outside of the energy spectrum of the 
system, we counted how many times the sign of the wave function changes over the length 
of the system for different values of E. Each new node corresponds to a new state of the 
system.^ 

For the random superlattice model we find it more convenient to use the phase formalism 
described, for instance, in Ref. 0. Within this formalism the density of state is expressed in 
terms of the phase variable, = tan -1 (ip'/ip). In the case of systems with a single band 
spectrum, this phase changes between and it when E sweeps the spectrum of the system 
from one band boundary to the other. In the superlattice, the spectrum of the wave in the 
absence of disorder consists of multiple bands. In this case, the phase increases by it across 
every allowed band, and stays constant and equal to nir, inside any n-th forbidden band. If 
disorder in our model is not too strong, the regions of the constant phase are preserved even in 
the presence of random fluctuations Fig. (fl|). and can be used for identifying the fluctuation 
boundaries of the bands in the disordered system. Then we can introduce a density of states 
N(E) for a single band, which is normalized to change from to 1, when energy spans from 
one fluctuation boundary to another. N(E) normalized this way is substituted in Eq. (JHJ) in 
order to calculate l s for the superlattice model. When disorder becomes stronger the regions 
of constant phase disappear, and the notion of the single band density of states becomes 
meaningless. In our calculations we always make sure to avoid such situations. 
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1 - 4 1 - 45 Wavenumber, k- a 1 " 55 1 - 6 

FIG. 1: The dependence of the phase near the band gap region (1.44 < ka < 1.56) separating the 
first and the second bands in the superlattice model, efo was taken from a uniform distribution 
with 5=0.1, L/a ~ 10 6 . 

III. SCALING DESCRIPTION OF THE MOMENTS OF THE DISTRIBUTION 
FUNCTION 



It was shown in Refs. 
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12 that the variance, a 2 of the Lyapunov exponent in the LLoyd 
model can be conveniently described in terms of a relationship between two scaling variables, 
t defined as 

r-^ (13) 

7 

and k, denned as 

k = ^ (14) 



In this paper we show that the variance of LE in more generic models can also be described 
in terms of the scaling function t(k). 

In order to demonstrate this result we computed a 2 and l s for different values of the energy, 
strength of disorder, and length of the system for both models under considerations. The 
results of these calculations were presented in the form of the function t(k), which is shown 
in Figs. 121 and 01 for the Anderson model and the superlattice model, respectively. The data 
included in these figures correspond to systems with L ^> l s , l\ oc . The first important result 
revealed by this figures is that all the data lie on a single curve, when expressed in terms 
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FIG. 2: Typical dependence of the scaling parameter r on k for the Anderson model. The width of 
the distribution of disorder changes from U = 0.08 to U = 0.16. Curves corresponding to different 
values of the width are not distinguishable. In the inset the region of small k is shown in the log-log 
scale. 
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FIG. 3: Dependence of the scaling parameter r on k for the dichotomic distribution (superlattice 
model) of d.2 with 5 = 0.1, 0.125, 0.15, 0.175, and 0.2. For every value of disorder we took 17 
lengths L, ranging from 320 to 20000 layers. Different offset value of Tu m was compensated (see 
text). In the inset the same is shown in the log- log scale. 
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of the variables r and k for both models. This result confirms our general conjecture that 
the second moment of the distribution function of LE can be universally described in terms 
of variables r and k regardless the microscopical nature of the models under consideration. 
While the form of the function r (/«), may differ for different models, its essential qualitative 
properties show a degree of universality: r (k) = 1 for k > 1, and it steeply decreases for 
k<C 1. We are most interested here in the latter region, where the fluctuation states arise. 
For the Lloyd model r = (tv/2) k for k <C l , 11 ' 12 while in the models studied in this paper 
the dependence of r upon k is much steeper. In order to obtain a better insight into the 
properties of r (k) for small k, we re-plotted our numerical data in log-log coordinates (see 
insets in Figs. |2]and|21). Before interpreting these figures we have to note that unlike the 
case of the Lloyd model, where r(0) = 0, in the models considered here t(0), while very 
small, is not equal to zero. The reason for this is the small fluctuations of the LE due to non- 
resonance tunnelling through a random barrier, which contributes to r at the fluctuation 
spectrum boundary where k = 0. This small contribution is model specific, and in the 
Anderson model it can be neglected everywhere with exception of a small neighborhood of 
the fluctuation spectrum boundaries. This can be seen from the fact that while k changes 
over at least two orders of magnitude, the data for the Anderson model (inset in Fig. |2J) 
form a straight line with exceptions of a few points corresponding to very small values of k. 
According to these results, t(k) has the form 

T=C K a +T lim , (15) 

where Tu m stands for the non-universal correction discussed above. In the superlattice model 
the value of ru m is more significant, and therefore has to be compensated. In order to 
estimate coefficients C and k, we select only those data for which l s < L and use linear 
regression. The results of the fit are presented in the table below: 





Anderson model 


Superlattice models 


c 


1.27 


1.08 


a 


0.27 


0.40 



These results demonstrate that while the nature of the scaling parameters is universal for 
both models, the numerical values of the respective parameters are model dependent. An 
interesting question is whether the values of C and a depend upon the type of statistics of 



11 



0.50 
0.45 
0.40 
0.35 
0.30 
8 0.25 
0.20 
0.15 
0.10 
0.05 



I ■ 



1.4 
1.3 
1.2 
1.1 
1.0 

0.9 U 

0.8 

0.7 

0.6 

0.5 



2 

LI I 



FIG. 4: Dependence of the index of the scaling parameter a (filled squares, left axis) and the factor 
C (circles, right axis) on L/l s for the Anderson model. 

the respective random parameters of our models (site energy for Anderson model, the layer 
width for the superlattice model). In the case of a superlattice model we found that the 
change in statistics (from the box to dichotomic distribution) did not affect the values of the 
coefficients C and a. For the Anderson model with the dichotomic distribution of the site 
energies the results were inconclusive. Strong noise in the data for the dichotomic process 
prevented us from positively establishing equivalency of the coefficients for the two different 
types of statistics. 

In the region of fluctuation states, a new intermediate regime of lengths L, in which 
koc <C L <C l s appears. This regime does not exist for in-band states. It is natural to 
anticipate that the scaling behavior of our systems in this regime would change. In order to 
study this question, we divided our data in groups according to the value of L/l s , including 
points with L/l s > 1 as well as with L/l s < 1. Carrying out statistical analysis of the data 
for fixed values of L/l s we were able to obtain dependencies of the parameters C and a on 
L/l s ; the respective results are presented in Figs. H] and 03 First of all, we would like to 
note that these dependencies saturate to values presented in the Table 1 for L/l s > 1. This 
confirms our assumption that in this regime r depends upon a single parameter, k. 

For shorter systems, however, a new parameter, L/l s emerges. For the Anderson model 
we were able to show that a(L/l s ) is best described by the logarithm a(L/l s ) ~ \n(l s /L), 
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10 20 30 40 50 

FIG. 5: The exponent of the scaling parameter a (filled squares, left axis) and the factor C (circles, 
right axis) as functions of L/l s for the superlattice model. Large and small symbols correspond to 
dichotomic and box distribution of d,2, respectively. 

which means that the variance of the Lyapunov exponent, a 2 , in this regime demonstrates 
an anomalous scaling with the length of the system, L: 

a 2 oc — — exp[a(L/Z s )lnK] oc L^ 1+lnK \ (16) 

Llloc 

It is interesting to note that when k decreases, 1 + In k may become negative resulting in 
a 2 increasing with L. This behavior can be qualitatively understood from the following 
arguments: The condition L <C l s means that for the most of the realizations of the random 
potential no states exist in the energy interval under discussion. The transmission through 
such realizations fluctuates rather weakly. The greatest contribution to the transmission 
fluctuations is given by those few realizations that can support at least a single state. The 
probability for such realizations to arise grows when the length of the system increases, 
resulting in the respective increase of a 2 . This behavior, of course, breaks down for very 
large values of l s , which correspond to states close to the genuine spectral boundary, because 
for these states a 2 is determined by a non- universal correction to r given by r^ m . 

The behavior of a 2 given by Eq. (J16|) can be confirmed by plotting directly the function 
<J 2 (L) for energies from the band-gap. Fig. |H1 presents such a plot for the Anderson model 
for the value of k equal to k — 0.2. It demonstrates a good agreement with Eq. (fTrJj) : the 
slope of the curve was found to be equal to 1.77, while an estimate for this slope from 
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-1.6 -1.4 -1.2 -1.0 -0.8 

FIG. 6: The logarithm of the scaling parameter t for the Anderson model as a function of log 10 L/l s 
for intermediate values of energy when l s is not too large. Points are the result of numerical 
calculations and the straight line is a linear fit. 

Eq. (fTljj) gives 1.78. It should be noted, however, that the regime described by Eq. (fTrJJ) 
exists in a relatively narrow interval of energies, at least for the Anderson model with the 
box distribution. The reason for this is that l s grows very fast in the region of fluctuation 
states when the energy is shifted toward the fluctuation spectrum boundary. Very large l s 
means that only few realizations of our system support at least a single state. Therefore, for 
the most realizations transmission occurs via non-resonant under-barrier tunnelling. The 
statistics of the transmission for this subset of realizations is determined by the localization 
length alone (l s is exact infinity for these realizations). As a result, we have a competition 
between a small number of realizations, supporting states, for which fluctuations of the 
Lyapunov exponent are large and grow with the length, and the majority of realizations, 
in which a 2 is small, and decreasing with length. At very large l s the contribution to a 2 
from the representative realizations becomes larger than the contribution from the resonant 
realizations, and Eq. (|16|) fails. In this case, an asymptotic behavior of a 2 is again controlled 
by the localization length alone, as it can be seen in Fig. where a 2 L saturates at L much 
smaller than l s . 

The assumption about the Gaussian form of the distribution of LE is the result of the 
central limit theorem, and strictly speaking is true only asymptotically when L — > oo. At 
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FIG. 7: The logarithm of the scaling parameter t for the Anderson model as a function of log 10 L/l s 
for energies corresponding to extremely large values of l s . The saturation occurs at the length close 
to the localization length. 



finite L the distribution function deviates from t 
SPS holdsji^^ However, it was found in Refs. 14 



le Gaussian form even in the regime when 
231 that this deviation, as measured by the 



magnitude of the third and higher moments, increases significantly in the vicinity of the 
band boundary of the initial spectrum. This result was obtained analytically for the white- 
noise potential in Ref. The first study of the scaling properties of the third moment was 



23 to the 



reported in Ref. |23j. In this part of the paper we expand scaling analysis of Ref. 
superlattice model, and compare the results obtained for these two models. We consider the 
scaling properties of the third cumulant g = ((7 — (7)) 3 ), which characterizes the asymmetry 
or skewness of the distribution function. Fig. |H] shows the energy dependence of the third 
moment for the Anderson model. It is seen that this moment significantly grows in the 
vicinity of the initial band boundaries of both models, which means that the significant 
deviation of the distribution function of LE from the Gaussian form in the region, where 
traditional SPS violates is a universal phenomenon. 

To analyze scaling properties of the third cumulant we consider the dimensionless param- 
eter 

r 3 = gL 2 l loc (17) 
The dependence of T3 on k for the Anderson model and superlattice is shown in Figs. 01 and 
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FIG. 8: Dependence of the renormalized third cumulant, gL 2 , on energy in the vicinity of the band 
edge of a pure system (U = 0.05) for the Anderson model. 
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FIG. 9: Dependence of the parameter t% = gL 2 li oc on k^ 1 = l s /li oc (Anderson model) for a set 
of different widths of the distribution of the potential: 0.001 < U < 0.21. Error bars show the 
dispersion of the results of numerical simulations near a mean values shown by squares. 

ITUl respectively. 

One can see from these figures that while data for the parameter r 3 are rather noisy, it 
shows a relatively good scaling behavior as a function of the single parameter k for both 
models. This fact itself is quite remarkable since it demonstrates that even in the region, 
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FIG. 10: Dependence of t 3 on li oc /l s for dichotomic distribution (superlattice model) of d,2 with 5 
= 0.1, 0.125, 0.15, 0.175, and 0.2. For every value of disorder we took 17 length L, ranging from 
320 to 20000 layers. On the insert the same is shown in log-log scale. 

where the distribution function of LE deviates significantly from the Gaussian form, it can 
still be characterized by two parameters within the scaling procedure suggested here. 

The better data quality for the superlattice model allowed for more thorough study of 
the third moment. The insert in Fig. 1101 shows a good scaling behavior similar to Eq. (|15|) : 

-T 3 = C 3 K aa +T 3>lim (18) 

The limiting value T 3 ji m , was substantially smaller than Tu m , so no explicit correction was 
needed to obtain Fig. [TUJ 

For intermediate lengths, l\ oc <C L <C l s , we analyzed data using approach similar to that 
employed to obtain Figs. E]and|3J For fixed values of L/l s we obtained dependencies of the 
parameters C 3 and a 2 on L/l s (Fig. (IT)) , and found the saturated values of C 3 = 0.73 and 
a 3 = 0.52 - the same for both dichotomic and box distributions. 

IV. COMPARISON WITH THE GAUSSIAN WHITE NOISE MODEL 

It is well known that under certain circumstances statistical properties of one-dimensional 
disordered systems in the vicinity of the band edges of the initial spectrum can be universally 
described by replacing an actual random potential by a Gaussian white noise potential.- 5 One 
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FIG. 11: 03 (filled squares, left axis), the exponent of T3, and factor C2 (circles, right axis) as 
functions of L/l s for superlattice model. Large and small symbols correspond to dichotomic and 
box distribution of d<i respectively. 



of the manifestations of this fact is that the statistical properties of LE in the Anderson model 
with the box distribution of the site energie a 14ll9l2 ° are very similar to those of the continuous 
model with the Gaussian white-noise potential, 5 and are characterized by the same scaling 
parameter E/D 2 ^ 3 , where E is the energy counted from the initial band boundary, and D 
is the variance of the random potential. It was noted in Ref. Q that the scaling parameter 
k is a single- valued function of the Gaussian scaling parameter E/D 2 / 3 for the white-noise 
model, so that in this case these two parameters are equivalent to each other. An important 
question now arises: whether the apparent universality of the scaling description, suggested 
in this paper, is a mere consequence of the fact that in the region of the fluctuation states all 
models can be reduced to the Gaussian model, or this universality reveals more fundamental 
properties of this spectral region. This question was partially discussed in the Ref. [l?], in 
which it was shown that the behavior of the second moment of the LE in the vicinity of E = 
of the Anderson model obeys the scaling description in terms of the parameter k, while the 
Gaussian approximation certainly does not work in this part of the spectrum. In this paper, 
we address this question considering regions of the fluctuation states in the superlattice 
model.— The inset to Fig. ^] shows the plot of the parameter r versus the Gaussian scaling 
parameter (k — ki)/D 2 ^ 3 (D oc 5 2 ), where ki is the dimensionless frequency of one of the 
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FIG. 12: Normalized variance of LE, r, plotted versus parameter {k — kj)/5 demonstrates a good 
scaling. Scaling with (k — fcj)/<5 4 / 3 , predicted by the Gaussian white noise model, shown in the 
inset, fails. The data was generated in the superlattice model with the box distribution for five 
values of 5 =0.1, 0.125, 0.15, 0.175, and 0.2, L=20000. We included the frequencies from two band 
edges - the upper edge of the first band and the lower edge of the second band. Altogether, the 
band gap region between the first and the second bands is covered entirely. 

initial band boundaries of the superlattice for several values of the disorder. Moreover, we 
included the frequencies from the upper edge of the first band and the lower edge of the 
second band. We found that instead of E/D 2 ^ 3 predicted by the Gaussian white noise model, 
our data are better scaled with the parameter E/D 1 ^ 2 . One can see from Fig. that while 
the Gaussian scaling fails, the function t(k) discussed in the previous section of the paper 
gives the best scaling description of this model as well as of the Anderson model. We can 
conclude, therefore, that the scaling parameter n retains its universal significance beyond 
the validity of the white-noise approximation. 

V. CONCLUSION 

In this paper we studied scaling properties of the distribution function of the Lyapunov 
exponent for two one-dimensional disordered models: the Anderson model with diagonal 
disorder, and the model of a scalar wave propagating in a random superlattice. The main 
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result of the paper is that in the region of band-edge and fluctuation states, where simple 
SPS fails, the distribution function can be described by two independent scaling parameters: 
the localization length, li oc , and an additional length l s , introduced in Refs. Illlll2l which is 
related to the integral density of states. The fact of principal importance is that not only 
the second moment of the distribution is described by these two parameters, but so also is 
the third moment. This means that even though in the region of fluctuation states the form 
of the distribution function strongly deviates from the Gaussian, it still can be described 
within the suggested two-parameter scaling approach. 

Among other important results of the paper we would like to note the detailed study of the 
properties of the variance and the third moment of LE in the region of fluctuation states. 
We showed that both, the normalized variance and the third cumulant presented by the 
scaling functions r and T3, demonstrate a power law dependence upon the scaling parameter 
k. Parameters of this power law dependence were found to depend weakly upon the type of 
statistics used to characterize our random systems, but are different for the Anderson model 
and the superlattice model. When the length of the system becomes smaller than l s , we 
showed that the scaling behavior of a 2 deviates significantly from the central limit theorem 
behavior even when L remains much bigger than the localization length. 
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